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Abstract. We present new HST NICMOS observations 
of NGC 4945, a starburst galaxy hosting a highly obscured 
active nucleus that is one of the brightest extragalactic 
sources at 100 keV. The HST data are complemented with 
ground based [Fen] line and mid-IR observations. 

A lOOpc-scale starburst ring is detected in Paa, while 
H2 traces the walls of a super bubble opened by supernova- 
driven winds. The conically shaped cavity is particularly 
prominent in Paa equivalent width and in the Pacv/H2 
ratio. Continuum images are heavily affected by dust ex- 
tinction and the nucleus of the galaxy is located in a highly 
reddened region with an elongated, disk-like morphology. 
No manifestation of the active nucleus is found, neither 
a strong point source nor dilution in CO stellar features, 
which are expected tracers of AGN activity. 

Even if no AGN traces are detected in the near-IR, 
with the currently available data it is still not possible to 
establish whether the bolometric luminosity of the object 
is powered by the AGN or by the starburst: we demon- 
strate that the two scenarios constitute equally viable al- 
ternatives. However, the absence of any signature other 
than in the hard X-rays implies that, in both scenarios, 
the AGN is non-standard: if it dominates, it must be ob- 
scured in all directions, conversely, if the starburst dom- 
inates, the AGN must lack UV photons with respect to 
X-rays. 

An important conclusion is that powerful AGNs can be 
hidden even at mid-infrared wavelengths and, therefore, 
the nature of luminous dusty galaxies cannot be always 
characterized by long-wavelength data alone but must be 
complemented with sensitive hard X-ray observations. 
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* Based on observations made with the NASA/ESA Hubble 
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1. Introduction 

A key problem in studies of objects emitting most of their 
energy in the FIR/submm is to establish the relative im- 
portance of highly obscured Active Galactic Nuclei (AGN) 
and starburst activity. In particular, it is important to 
know if it is still possible to hide an AGN, contributing 
significantly to the bolometric emission, when optical to 
mid-IR spectroscopy and imaging reveal only a starburst 
component. 

Several pieces of evidence suggest that most cosmic 
AGN activity is obscured. Most, and possibly all, cores of 
large galaxies host a supe rmassive black hole (10 6 -10 9 M Q ; 
e.g. Richstone et al. 1998| ). To complete the formation pro- 
cess in a Hubble time, accretion must proceed at high 
rates, producing quasar luminosities (L ~ 10 12 L Q ). How- 
ever the observed black hole density is an order of magni- 
tude greater than that expected from the observed quasar 
light, assuming accretion efficiency of 10%, suggesting that 
most of the ac cretion history is obscured (e.g. Fabian & 
Iwasawa 1999, and references therein). It is estimated ei- 
ther that 85% of all AGNs are obscured (type 2) or that 
85% of the accretion history of an object is hidden from 
view. 

In addition, the hard X-ray background (> 1 keV) re- 
quires a large population of obscured AGNs at higher red- 
shifts (z ~ 1) since the observed spectral energy distribu- 
tion cannot be explained with the continua of Quasars, i.e. 



un-obscured (type 1) AGNs (Comastri et al. 1995, Gilli 
et al. 1999). Despite the above evidence, detections of ob- 
scured AGNs at cosmological distances are still sparse (e.g. 
Akiyama et al. |1999[ ). 

Ultra Luminous Infrared Galaxies (ULIRGs; 



see 



Sanders & Mirabel 1996 for a review) and the sources 
detected in recent far-infrared and subrrrm surveys per- 
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formed with ISO and SCUBA 
al 



(e.g. Rowan-Robinson et Table 1. Log of HST observations. 



1997, Blain et al. 199E and references therein) are can- 
didate to host the missing population of type 2 AGNs. 
However, mid-IR ISO spectroscopy has recently shown 
that ULIRGs are mostly powered by starbursts and that 
no trace of AGNs is found in the majority of cases (Gen- 
zel et al. 1998| ; Lutz et al. 1998| ). Yet, the emission of a 
hidden AGN could be heavily absorbed even in the mid- 
IR. Indeed, the obscuration of the AGN could be related 
to the starburst phenomenon, as ob served for Seyfert 2s 
(Maiolino et al. |l995| ). Fabian et al. (|l998|) proposed that 
the energy input from supernovae and stellar winds pre- 
vents interstellar clouds from collapsing into a thin disk, 
thus maintaining them in orbits that intercept the major- 
ity of the lines of sight from an active nucleus. 

In this paper, we investigate the existence of com- 
pletely obscured AGNs and the Starburst- AGN connec- 
tion through observations of NGC 4945, one of the closest 
galaxies where an AGN and starburst coexist. NGC 4945 
is an edge-on (i ~ 80°), nearby (D — 3.7Mpc) SB spiral 
galax y hosting a powerful nucle ar starburst (Koornneef 
1993|; Moorwood & Oliva 1994a ). It is a member of the 



Centaurus group and, like the more famous Centaurus A 
(NGC 5128), its optical image is marked by dust extinc- 
tion in the nuclear regions. The ONLY evidence for a hid- 
den AGN comes from the hard X-rays where NGC 4945 is 
characterized by a Compton-thick spectrum (with an ab- 
sorbing column density of A^h = 5 x 10 24 cm~ 2 , Iwasawa 



et al. 1993) and one of the brightest 100 keV emissions 



among extragalactic sources (Done et al. 1996). Recently, 
BeppoSAX clearly detected variability in the 13-200 keV 
band (Guainazzi et al., 2000). 

Its total infrared luminosity deriv ed from IRAS data 
is ~ 2.4 x 10 10 L Q (Rice et al. |l988f ), ~ 75% of which 
arises from a regio n of < 12" x 9" centered on the nucleus 
(Brock et al. 1988 ). Although its star formation and super- 
nova rates are mode rate, ~ O.4M yr _1 and ~ 0.05 yr -1 



(Moorwood & Oliva ^994aj), the starburst activity is con- 
centrated in the central ~ 100 pc and has spectacular con- 
sequences on the circumnuclear region which is character- 
ized by a conical cavity evacuated by a supernova-driven 



wind (Moorwood et al. 1996a). 

The radio emission is characterized by a compact non- 
thermal core with a luminosity of ~ 8 x 10 38 ergs _1 (El- 
mouttie et al. 1997 ). It is one of the first H2O an d OH 
mega maser sources detected (dos Santos & Lepinc 1979| ; 
Baan 1985[) and the H2O maser was mapped by Greenhill 
et al. ( 1997 ) who found the emission linearly distributed 
along the position angle of the galactic disk and with a ve- 
locity pattern suggesting t he pr esence of a ~ 10 6 M Q black 
hole. Mauersberger et al. ( |1996| ) mapped the J = 3 — 2 line 
of 12 CO which is mostly concentrated within the nuclear 
~ 200 pc. 

We present new line and continuum images obtained 
with the Near Infrared Camera and Multi Object Spectro- 
graph (NICMOS) on-board the Hubble Space Telescope 



Dataset 


Filter 


Texp (s6c) 


Description 


n4mq01010 


F110W 


768 


.1 


n4mq02010 


F160W 


768 


H 


n4mq01040 


F222M 


288 


K 


n4mq01070 


F237M 


288 


CO 


n4mqal010 


F222M 


288 


K background 


n4mqaiU2U 


r Z61NL 


288 


CO background 


n4mqblnrq 


F187N 


320 


Paa 


n4mqblnuq 


F190N 


320 


Paa continuum 


n4mqblo0q 


F190N 


320 


Paa continuum 


n4mqblo3q 


F187N 


320 


Paa 


n4mqblo9q 


F187N 


320 


Paa 


n4mqblodq 


F190N 


320 


Paa continuum 


n4mqbloiq 


F190N 


320 


Paa continuum 


n4mqblolq 


r 10 (IN 


o2U 


Paa 


n4mqblopq 


F212N 


320 


H 2 


n4mqblosq 


F190N 


320 


H2 continuum 


n4mqblovq 


F190N 


320 


H2 continuum 


n4mqblozq 


F212N 


320 


H 2 


n4mqblp2q 


F212N 


320 


H 2 


n4mqblp5q 


F190N 


320 


H 2 continuum 


n4mqblp9q 


F190N 


320 


H 2 continuum 


n4mqblpcq 


F212N 


320 


H 2 


u29r2p01t 


F606W 


80 


R archive 


u29r2p02t 


F606W 


80 


R archive 


u2e67z01t 


F606W 


500 


R archive 



(HST), aimed at detecting AGN activity in the near- 
infrared. These observations are complemented by recent 
ground based near- and mid-IR observations obtained at 
the European Southern Observatory. Section || describes 
the observations and data reduction techniques. Results 
are presented in Section || and discussed in Section |[ Fi- 
nally, conclusions will be drawn in Sec. ^. Throughout the 
pa per w e assume a distance of 3.7 Mpc (Mauersberger et 
al. 1996), whence 1" corresponds to ~ 18 pc. 



2. Observations and Data Reduction 

The nuclear region of NGC 4945 was observed on March 
17 th and 2 5 th , 1998, with NICMOS Camera 2 (MacKenty 
et al. |l997t ) using narrow and broad band filters for imag- 
ing in lines and continuum. HST observations are logged 
in Tab. [I]. All observations were carried out with a MUL- 
TIACCUM sequence (MacKenty et al. |1997| ) and the de- 
tector was read out non-destructively several times dur- 
ing each integration to remove cosmic rays hits and cor- 
rect saturated pixels. For each filter we obtained several 
exposures with the object shifted by ~ 1" on the detec- 
tor to remove bad pixels. The observations in the F222M 
and F237M filters were also repeated on a blank sky area 



Marconi et al.: The elusive AGN of NGC 4945 



3 





Fie. 2. Ca.1 Pan- imaffn. Svmhols a.s in Fie\ 1. The blank contours are from the Ho<-I-[Ntt1 imace hv Moorwoorl et a.1. 
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Table 2. Log of Ground Based observations. 



Image T exp (min) Date Instrument Telescope 



L' 10 May 30, 1996 IRAC1 ESO/MPI 2.2m 

N 40 May 27, 1996 TIMMI ESO 3.6m 

[Fen] 24 April 1, 1998 IRAC2B ESO/MPI 2.2m 



several arcminutes away from the source to remove ther- 
mal background emission. For narrow band images, we ob- 
tained subsequent exposures in line and near continuum 
filters with the object at several positions on the detector. 

The data were re-calibrated u sing t he pipeline software 
CALNICA v3.2 (Bushouse et al. |l997|) . A small (few per- 
cent) drift in the NICMOS bias level caused an error in the 
flat-fielding procedure which resulted in spurious artifacts 
in the final images (the s o-call ed "pedestal problem" - 
Skinner, Bergeron & Daou 1998 ). Given the strong signal 
from the galaxy, such artifacts are only visible in ratio or 
difference images. This effect was effectively removed using 
the pedestal estimation and quadrant equalization soft- 
ware developed by Roeland P. van der Marel which sub- 
tracts a constant bias level times the flat-field, minimizing 
the standard deviation in the images. For each filter, the 
corrected images were then aligned via cross-correlation 
and combined. Flux calibration of the images was achieved 
by multiplying the count rates (adus -1 ) for the PHOT- 
FLAM ( erg cm -2 adu -1 ) conversion factors (MacKenty et 
al. [1997| ). 

The narrow band images obtained at wavelengths ad- 
jacent to the Paa and H2 lines where used for continuum 
subtraction. The procedure was verified by rescaling the 
continuum by up to ±10% before subtraction and estab- 
lishing that this did not significantly affect the observed 
emission-line structure. 

WFPC2 observations in the F606W (R band) filter 
were retrieved from the Hubble Data Archive and 
re-calibrated with the standard pipeline software (Biretta 
et al. |1996[ ). 

Ground-based observations were obtained at the Euro- 
pean Southern Observatory at La Silla (Chile) in the con- 
tinuum L' (3.8 fim) and N (10 /im) bands, and in the [Fe 11] 
1.64 /j,m emission line and are logged in Tab. |2|. The L' im- 
age was obtained with IRAC1 (Moorwood et al. 1994b| ) at 
the ESO/MPI 2.2 m telescope on May 30, 1996 using an 
SBRC 58x62 pixel InSb array with a pixel size of 0'.'45. 
Double beam-switching was used, chopping the telescope 
secondary mirror every 0.24 s and nodding the telescope 
every 24 s to build up a total on-source integration time 
of 10 minutes in a seeing of 0'.'9. The N-band image was 



[Fen] 1.64 /xm image was taken with with the IRAC2B 



obtained with TIMMI (Kaufl et al. 1994) at the ESO 3.6m 



camera (Moorwood et al. 1992) on the ESO/MPI 2.2m 



telescope on April 1, 1998, using a 256x256 Rockwell NIC- 
MOSS HgCdTe array with 0'.'51 pixels. The [Fen] line was 
scanned with a A/AA = 1500 Fabry-Perot etalon cover- 
ing three independent wavelength settings on the line and 
two on the continuum on either side of the line, for a total 
integration times of 24 minutes on the line in 0'.'9 seeing. 
Standard procedures were used for sky subtraction, flat 
fielding, interpolation of hot and cold pixels at fixed posi- 
tions on the array, recentering and averaging of the data. 
For the [Fen] data, the continuum was determined from 
the two off-line channels and subtracted from the on-line 
data. The integrated [Fe 11] line flux is in excellent agree- 
ment with the value determined by Moorwood and Oliva 



(1994a) 



3. Results 

3.1. Morphology 

Panels a-d in Figure [j] are the continuum images in the 
NICMOS K, H, J and WFPC2 R filters^ The cross marks 
the position of the K band peak and the circle is the posi- 



tion of the H2O maser measured by Greenhill et al. (1997). 
The radius of the circle is the ±1" r.m.s. uncertainty of 
the astrometry performed on the i mages and based on the 
Guide Star Catalogue (Voit et al. 1997 ). The position of 
the K peak is offset by ~ 0'.'5 from the location of the 
H2O maser, hereafter identified with the location of the 
nucleus of the galaxy. Note that this offset is still within 
the absolute astrometric uncertainties of the GSC and the 
K peak could be coincident with the nucleus. The con- 
tinuum images are also shown with a "true color" RGB 
representation in Fig. |f (Red=F222M, Green=F110W, 
Blue=F606W). A comparison of photometry between our 
data and earlier published results is not straightforward 
since the NICMOS filters are different from the ones com- 
monly used. However, as shown in table ||, our measured 
fluxes in 6" and 18" circular apertures centered on the K 
band peak are within 15-30% of the ones by Moorwood & 



Glass (1984) measured in the same areas. 

Figure |l e is the H-K color map. H-K contours are also 
overlayed on the J image in Fig. [j]c. At the location of 
the maser, South-East of the K band peak, emission from 



telescope on May 27, 1996 using a 64x64 Si:Ga array with 
0'.'46 pixels. Again using double beam switching, total on- 
source integration time was 40 minutes in 1" seeing. The http://www.arcetri.astro.it/~marconi 
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Fig. 3. (a) L band contours overlayed on the NICMOS K band image (displayed with a logarithmic look-up table). 
Contours are 0.005 and from 0.01 to 0.14 with step of 0.01 (units of 10~ 16 erg cm~ 2 s _1 A _1 ). The frame boxes are 
centered on the nucleus, identified with the H2 maser position, (b) N band contours overlayed on the NICMOS K 
band image. Contours are from 0.02 to 0.08 with step of 0.003 (units and notation as above). 



Table 3. 

data. 



Comparison with ground based photometric 



Band 


(") 


Moorwood 
mag 


& Glass 1984 
Flux" 


This work 
Flux" 


K 


6" 


9.34 


7.1xl0 -15 


9.5xl0" 15 


K 


18" 


8.12 


2.2xl0~ 14 


2.5xl0~ 14 


H 


6" 


10.7 


5.7xl0~ 15 


7.1xl0~ 15 


H 


18" 


9.15 


2.4xl0~ 14 


2.7xl0~ 14 


J 


6" 


12.70 


2.4xl0~ 15 


2.7xl0~ 15 


J 


18" 


10.80 


1.4xl0 -14 


1.6xl0~ 14 



a In units of ergs cm A. 

The K band of ground based observations corresponds to 
the F222M NICMOS filter. Similarly H and J corresponds to 
F160W and F110W, respectively. 



galactic stars is obscured by a dust lane oriented along 
the major axis of the galactic disk. The morphology of 
this resembles an edge-on disk with a 4'.'5 radius (80 pc) 
and probably marks the region where high density molec- 
ular material detected in CO is concentrated (e.g. Mauers- 
berger et al. 1996). The average H-K color is ~ 1.7 with 
The dust lane has a sharp southern 



a peak value of 2.3. 
edge which is not very evident in the color map and can be 
explained by saturated absorption: background H and K 



emission becomes completely undetectable and the color 
is dominated by foreground stars. Therefore, the sharp K 
edge is evidence for a region with such high extinction 
that it is not detected even in the near-IR. With this dust 
distribution, the observed morphology in the continuum 
images is the result of an extinction gradient in the direc- 
tion perpendicular to the galactic disk. Patchy extinction 
is also present all over the field of view. At shorter wave- 
lengths, the morphology is more irregular because dust 
extinction is more effective (the same effect seen so ob- 
viously in Centaurus A, cf. Schreier et al. 199S, Marconi 



et al. 2000) and the conical cavity extensively mapped by 
Moorwood et al. (1996a) becomes more prominent: there, 
the dust has been swept away by supernova-driven winds. 
Indeed, in the R image, significant emission is detected 
only in the wind-blown cavity which presents a clear con- 
ical morphology with well defined edges and apex lying 
~ 3" from the K peak. Due to the above mentioned red- 
dening gradient, the apex of the cone gets closer to the 
nucleus with increasing wavelength (compare with R band 
and Paa/H2 images - see below). 

The continuum-subtracted Paschen a image (Fig. ||a) 
shows the presence of several strong emission line knots 
along the galactic plane, very likely resulting from a cir- 
cumnuclear ring of star formation seen almost edge-on. 



"Knot B" of Moorwood et al. (1996a) is clearly observed 
South East of the nucleus while " Knot C" , North- West of 
the nucleus, is barely detected. Both knots are also marked 
on the figure. 
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Dust extinction strongly affects the Paa morphology 
making very difficult to trace the ring and locate its center; 
a likely consequence is the apparent misalignment between 
the galaxy nucleus and the ring center. The observed ring 
of star formation is similar to what has been found in other 



starburst galaxies (cf. Moorwood 1996b). The starburst 
ring could result from two alternative scenarios: either the 
starburst originates at the nucleus, and then propagates 
outward forming a ring in the galactic disk; or the ring 
corresponds to the position of the inner Lindblad reso- 
nance where the gas density is naturally incre ased by flow 
from both sides (see the review in Moorwood 19961: ). 

Panel b shows the continuum-subtracted H2 image 
which traces the edges of the wind-blown cavity. As ex- 
pected, the morphology is completely different from that 
of Paa which traces mainly starburst activity. Note the 
strong H 2 emission close to the nucleus at the apex of 
the cavity with an elongated, arc-like morphology. The H2 
flux in a 6" x 6" aperture centered on the K band peak is 
1.1 x 1CP 13 erg cm -2 s _1 and corresponds to ~ 70% of the 
total integrated emission in the NICMOS field of view. 
This is in good agre ement with the 1.29±0.05 found by 
Koornneef & Israel ( 1996; ) in an equally sized aperture 



and the integrated 3.1x10 13 erg cm 2 s 1 from the map 



by Moorwood & Oliva ( 1994a ). We remark that contam- 
ination of H2 emission by HeiA2.112/im is unlikely since 
the line was detected neither by Koornneef (1993) nor by 
Moorwood & Oliva (1994a) and from their spectra we can 
set an upper limit of 5-10% to the Hei/H2 ratio. 

Panel c in Figure || shows that the equivalent width 
of Paa is up to 150-200 A in the star forming regions, 
but much lower in the wind-blown cavity. Since near-IR 
continuum emission within the cone is not significantly 
higher than in the surrounding medium, the low equivalent 
width within the cone is due to weaker Paa emission, the 
likely consequence of low gas density. 

Panel d in Figure || is the Paa/H2 ratio image which 
also traces the wind-blown cavity. Note that the cone 
traced by Paa and H2 is offset with respect to the light 
cone observed in R: this is a result of the reddening gra- 
dient in the direction perpendicular to the galactic plane. 

Fig. ||f is a true color RGB representation of line and 
continuum images (Red=F222M, Green=H 2 , Blue=Paa). 

L' and N band ground based images are shown in Fig. 
H a and b with contours overlayed on the NICMOS K band 
image. No obvious point source is detected at the location 
of the nucleus and the extended emission is smooth and 
regular, elongated as the galactic disk. 

The [Fe 11] emission shown in contours in Fig. ^|b devi- 
ates from the Paa image in a number of interesting ways. 
First, the northern edge of the cavity outlined most clearly 
in H2 emission is also detected, although more faintly, in 
[Fen], presumably excited by the shocks resulting from 
the superwind. Otherwise, the [Fen] emission displays two 
prominent peaks in the starburst region traced by Paa, 
one peak close to the nucleus and one offset at a position 



angle of about 250 degrees (counterclockwise from North) . 
In both of these regions the [Fen] /Paa ratio is much 
higher than in the rest of the starburst region. The [Fe 11] 
emission likely originates in radiative supernova remnants 
(SNRs). In the dense nuclear region of NGC4945 the ra- 
diative phase of the SNRs will be short, and hence the 
[Fe 11] emission will be much more strongly affected by the 
stochastic nature of supernova explosions in the starburst 
ring than Paa. The regions of high [Fen] /Paa ratios thus 
simply trace recent supernova activity. 

3.2. Reddening 

A lower limit and a reasonable estimate of reddening can 
be obtained from the H-K color image in the case of fore- 
ground screen extinction. In this case, the extinction is 
simply 



Ay = 



E(H - K) 
<H) - c(K) 



(1) 



where the color excess is given by the difference between 
observed and intrinsic colour, E(H — K) = (H — K) — 
(H — K) and the c coefficients represent the wavelength 
dependence of the extinction law; A\ — c(X)Ay. We 
have assumed A\ = A\ Mm (A/l jam) -1,75 (A > 1 /zm) and 
Ay = 2A2Ai lim . Spiral and elliptical galaxies have aver- 
age intrinsic colours (H — K) Q ~ 0.22 with 0.1 mag dis- 
persion (Hunt et al. 1997) and the color correction due to 
the non-standard filters used by NICMOS is negligible - 
(H — K)o ~ 0.26 instead of 0.22. In the region of the Paa 
ring, the average color H — K = 1.1 yields Ay ~ 11, in fair 
agreement with the estimate Ay > 13 from the Balmer 
decrement presented below. In knot B (H — K) = 1.2 
yields Ay = 12.5, while in knot C (H - K) = 0.64 yields 
Ay = 5.2. 

A different reddening estimate can be derived from 
the analysis of Hydrogen line ratios. We can estimate 
the reddening to "Knot B" and "Knot C" by using 
t he im ages and spectra published by Moorwood et al. 
(1996a). The inferred reddening (assuming an intrinsic 
ratio Paa/Ha=0.18, and A(Ra) = 0.81A V , A(Paa) = 
0.137A V ) is Ay=3.2 for Knot B and A v =3.8 for Knot C. 
We can also estimate a lower limit to the reddening on 
the Paa ring. Considering a region ~ 11" x 5" aligned 
along the galactic plane, including all the stronger Paa 
emission, we find Paa/Ha > 500 which corresponds to 
Ay > 13mag, a value in agreement with the estimate 



given by Moorwood & Oliva ( |1988j ), Ay = 14 ± 3, from 
the Bra/Br/7 ratio in a 6" x 6" aperture centered on the 
IR peak. 

We note that the first approach measures the mean 
extinction of the starlight, while the second one measures 
the extinction toward the HII regions. Therefore, these Ay 
estimates indicate that in the case of Knot C the star light 
and the emitting gas are located behind the same screen. 
Conversely, Knot B has a lower extinction and must there- 
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fore be located in front of the screen hiding the star light. 
A likely interpretation is that Knot C is located within 
the galactic plane on the walls of the cavity farthest from 
us. whereas Knot B is located above the galactic plane, 
toward the observer 

It appears that the hypothesis of screen extinction can 
provide reasonable results. Of course the true extinction, 
i.e. the optical depth at a given wavelength, is larger if dust 
is mixed with the emitting regions. However, it should 
be noted that the case in which dust is completely and 
uniformly mixed with the emitting regions does not apply 
here because the observed color excesses are larger than 
the maximum value expected in that case (E(H — K) ~ 
0.6). 

3.3. CO Index 

A straight computation of the CO stellar index as 
W{CO) = m(CO) - m(K), where m{CO) and m{K) are 
the magnitudes in the CO and K filters, is hindered by the 
high extinction gradients. Therefore we have corrected for 
the reddening using the prescription described above: 

W(CO) = m(CO) - m{K) + 2^_^1 E{H - K){2) 

where, as above, the c coefficients represent the wave- 
length dependence of the reddening law. The correction 
is 0.145 E(H — K) which is important since the expected 
CO index is ~ 0.2. 

The "corrected" photometric CO index map is dis- 
played in Fig. He. 

As a check, in the central 4" x 4" we derive a photo- 
metric CO index of 0.18 which is in good agreement with 
the value 0.22 obtained from spectroscopic observations 
by Oliva et al. ( 1995] ), when one takes into account the 
uncertainties of reddening correction. 

In the central region there are three knots where the 
CO index reaches values ~ 0.25 aligned along the galactic 
disk. However, we do not detect any clear indication of 
dilution by a spatially unresolved source, that would be 
expected in the case of emission by hot (~ 1000 K) dust 
heated by the AGN. There are regions close to the location 
of the H2O maser where the CO index is as low as 0.08 
but that value is still consistent with pure stellar emission 
or, more likely, with an imperfect reddening correction. 

3.4. AGN activity 

The NICMOS observations presented in this paper were 
aimed at detecting near-IR traces of AGN activity in the 
central (R < 10") region of NGC 4945. Indeed, recent 
NICMOS studies exploiting the high spatial resolution of 
HST show that active galactic nuclei are usually charac- 
terized by prominent point sources in K, dete cted e .g. in 
the Seyfert 2 galaxies Circinus (Maiolino et al. 1999a ) and 
NGC 1068, and the radio galaxy Centaurus A (Schreier et 



al. 1998). NGC 4945 does not show any point-like emis- 
sion at the position of the nucleus (identified by the H 2 
maser) and the upper limit to the nuclear emission is 
F A (F222M) < 2 x 10~ 13 erg cm" 2 s" 1 /im" 1 . 

We also do not detect any dilution of the CO absorp- 
tion features by hot dust emission, as observed in many 
active galaxies (Oliva et al. 1999b| ). From the analysis of 
the CO index image, non-stellar light contributes less than 
F A (F222M) < exlO-^ergcm^s-Vm^ 1 thus providing 
a tighter upper limit than above. 

The lack of a point source in the ground based L and N 
observations also places upper limits on the mid-IR emis- 
sion, though less tight due to the lower sensitivity and 
spatial resolution (F A (L) < 1.2 x 10 -12 ergcm~ 2 s _1 /1m" 1 
and F A (N) < 6.0 x 10~ 13 erg cm" 2 s - V m_1 )- 

Finally, type 2 AGNs are usually characterized by ion- 
ization cones detected either in line images or in exci- 
tation maps, i.e. ratios between high and low excitation 
lines (usually [Om] and Ha) revealing higher excitation 
than the surrounding medium. In NGC 4945 the equiv- 
alent width of Paa and the Paa/H2 ratio indeed show a 
cone morphology but the behaviour is the opposite of what 
expected, i.e. the excitation within the cone is lower than 
in the surroundings and the H2/Paa ratio increases up to 
~ 5 (see Fig. ||d). Two processes could be responsible for 
the enhanced H2 emission - either shocks caused by the 
interaction between the supernova-driven wind and the 
interstellar medium or exposure to a strong X-ray domi- 
nated photon flux emitted by the AGN. But in any case 
there is absolutely no indication of the strong UV flux 
which produces "standard" AGN ionization cones. 

We find, therefore, no evidence for the expected AGN 
markers in our NICMOS data. 



4. Discussion 

Although no trace of its presence has been found in these 
data, the existence of an obscured AGN in the nucleus of 
NGC 4945 is unquestionably indicated by the X-rays (Iwa- 
sawa et al. 1993, Done et al. 1996). Recent, high signal-to- 
noise observations by BeppoSAX (Guainazzi et al., |2000[ ) 



have confirmed the previous indications of variability from 



Ginga observations (Iwasawa et al. 1993): in the 13-200 



keV band, where the transmitted spectrum is observed, 
the light curve shows fluctuations with an extrapolated 
doubling/halving time scale of r ~ 3 — 5 x 10 4 s. These 
time scales and amplitudes essentially exclude any known 
process for producing the high energy X-rays other than 
accretion onto a supermassive black hole. 

Making the 3xl0 42 ergs _1 observed in the 2-10 keV 
band with BeppoSAX would require about 10000 of the 
most luminous X-ray binaries observed in our Galaxy (e.g. 
Scorpio-Xl) and only a few of this objects are known. 

Alternatively, very hot plasma (KT ~ a few keV) , due 
to supernovae, has been observed in the 2-10 keV spectrum 
of starburst galaxies, but at higher energies (> 30 keV) the 
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emission is essentially negligible (Cappi et al. 1999; Persic 
et al. 1998); whereas the emission of NGC4945 peaks be- 
tween 30 and 100 keV. Also, given that the X-ray emission 
is observed through a gaseous absorbing column density 
of a few times 10 24 cm -2 , both the 10000 superluminous 
X-ray binaries and the very hot SN wind should be hid- 
den by this huge gaseous column. It is very difficult to 
find a geometry for the gas distribution that could pro- 
duce this effect. We therefore conclude that the presence 
of an AGN provides the only plausible origin of the hard 
X-ray emission. 

The above considerations combined with the absence 
of any evidence for the presence of an AGN at other wave- 
lengths has important consequences irrespective of the rel- 
ative, and unknown, contributions of the starburst and 
AGN to the total bolomctric luminosity. This is illustrated 
below by considering the extreme possibilities that the lu- 
minosity is dominated either by the starburst or the AGN. 

4-1- NGC 4945 as a starburst dominated object 

Most previous studies have concluded that the FIR emis- 
sion in NGC 4945 can be attributed solely to starburst 



activity (e.g. Koornnccf 1993, Moorwood & Oliva 1994a) 
without invoking the presence of an AGN. 

We note that, on average, active galaxies are char- 
acterized by LfirI 'Lsr-y ratios much larger than star- 
bursts and this fact was sometimes invoked to discern 
starbursts from AGNs (see the discussion in Genzel et 
al. 1998). In this regard, NGC 4945 has a starburst-like 
ratio: I<FiR,/£Br7 ~ 1.4 X 10 5 (from observed Paa with 
j4v=15mag). This is similar to the value for the prototyp- 
ical starburst galaxy M82 (Lfir/Lbtj ~ 3.4 x 10 5 , Rieke 



et al. 198C| ), suggesting that the FIR emission of NGC 



4945 may arise from the starburst. 



Genzel et al. (1998) showed that, when considering the 
reddening correction derived from the mid-IR - usually 
much larger than from the optical and near-IR - the ob- 
served H line emission from the starburst translates into 
an ionizing luminosity comparable to the FIR luminosity. 
Indeed, if in NGC 4945 the bulk of H emission is hidden 
by just j4y=45mag, LFm/iion ~ 1 and the observed star- 
burst activity is entirely responsible for the FIR. 

Although all the bolometric luminosity could be gen- 
erated by a starburst it is also possible to construct star- 
burst models which are consistent with the observed near 
infrared properties but generate a much lower total lumi- 
nosity. It is important to recall that Lfir/I/b^ represents 
the ratio between star formation rates averaged over two 
different timescales, i.e. > 10 8 yrs and < 10 7 yrs, respec- 
tively. Therefore, this ratio strongly depends on the past 
star formation history. For example, objects which have 
not experienced star formation in the past 10 7 yrs will emit 
little Br7, but significant FIR radiation. A more quantita- 
tive approach is presented in Fig. ^ where we compare the 
observed nuclear properties of NGC 4945 with synthesis 



models by Leitherer et al. (1999). We have considered two 
extreme cases of star formation history. The thick solid 
line in the figure represents an instantaneous burst with 
mass 3.5 x 1O 7 M0 whereas the thick dashed line is a con- 
tinuous star formation rate of O.13M0 yr -1 . In both cases 
a Salpeter initial mass function (i.e. oc M~ 2 - 35 ), upper 
mass cutoff of 100M Q and abundances Z = Z Q are cho- 
sen. Panel 1 shows the evolution of the ionizing photon 
rate (Q(H)) as a function of time after the beginning of the 
burst. The shaded region limits the values compatible with 
the observations; Q(H) is estimated from the total Paa 
flux in the NICMOS images (5.6 x 10~ 13 erg s -1 cm" 2 ), 
dereddened with Ay = 5mag and Ay — 20mag and con- 
verted using case B approximation for H recombinations. 
Panel 2 gives the equivalent width of Bi'7 (Wx(Br7)); the 
observed value outlined by the shaded area is a lower limit 
for the starburst models and was derived by rescaling the 
observed Paa flux and dividing by the flux observed in the 
same aperture with the F222M filter. Panel 3 is the evolu- 
tion of the Super Nova Rate (SNR). Estimates of SNR from 
radio observations suggest values > 0.3 yr~ x (Koornneef 
|1993[ ), 0.2 yr" 1 (For bes fc N orris |1998| ), down to 0.05 yr" 1 
(Moorwood & Oliva 1994a). The shaded region covers the 



0.01-0.4 yr range. Panel 4 is the mechanical luminosity 
produced by the Supernovae. Finally, panels 5 and 6 give 
the K-band and bolomctric luminosity, respectively. The 
allowed range for the K monochromatic luminosity is given 
by the total observed flux in a 6" x 6" aperture centered on 
the K peak where photospheric emission from supergiants 



is known to dominate (Oliva et al. 1999b). The upper and 



lower limits represent the values obtained after deredden- 
ing by Ay = 5mag and Ay — 20mag. The upper limit to 
the bolometric luminosity is the total NGC 4945 luminos- 
ity derived from IRAS observations (Rice et al. 198S). In 
all cases the thin dotted line represents the time at which 
the properties of the instantaneous burst meet the ob- 
servational constraints. The crossed square represent the 
combination of the two models at t — 10 74 yr. 

It is clear from the figure that an instantaneous burst 
of t ~ 10 6,8 yr is capable of meeting all the observational 
constraints. It reproduces the correct supernova rate and 
K band luminosity and its bolometric luminosity domi- 
nates the total bolometric luminosity of the galaxy. Con- 
versely the continuous burst fails to reproduce the SNR 
and K luminosity. Just considering these two models alone 
it is tempting to infer that the starburst powers the bolo- 
mctric emission of NGC4945. However, the instantaneous 
and continuous SFR are two extreme and simplistic cases. 
More realistically the SF history is more complex since 
bursts have a finite and limited length or are the combi- 
nation of several different events. As an example we con- 
sider the case of two bursts of star formation taking place 
at the same time: one instantaneous and the other contin- 
uous. Both have the same characteristics as the bursts pre- 
sented above. The properties of this double burst model at 
t = 10 7,4 yr are shown in the figure by the crossed squares. 
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Fig. 4. Properties of a burst of star formation as a function of the time elapsed from the beginning of the burst 
(models by Leitherer et al. 1999). The thick solid line represents an instantaneous burst with mass 3.5 x 10 7 M Q . The 
thick dashed line is a continuous star formation rate of 0.13Mq yr -1 . See text for more details on the models. Panel 1 is 
the time dependence of the ionizing photon rate. The shaded area limits values consistent with observations. The thin 
dotted line is drawn at a time in which the "instantaneous" burst meets the observational constraints. The crossed 
square represents the combination of the properties of the two models at t = 10 7 4 yr. Panel 2 gives the Br7 equivalent 
width. The shading outline the lower limit given by observations. The other symbols are as in panel 1. Panel 3 gives 
the Supernova Rate. As above the shaded area limits the range of values allowed from observations. Symbols as in 
panel 1. Panel 4 gives the mechanical luminosity. Symbols as in panel 1. Panel 5 gives the monochromatic luminosity 
in the K band (ergs _1 A). Symbols as in panel 1. Panel 6 gives the bolometric luminosity of the burst. The shading 
marks the upper limit set by the total IRAS luminosity of the galaxy. 
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The choice of the time is arbitrary and any other value be- 
tween 10 7 ' 2 yr and 10 75 yr might do. Even in this case the 
starburst model meets all the observational constraints: 
Q(H) is provided for by the continuous burst while SNR 
and K luminosity come from the instantaneous burst. The 
important difference with respect to the single instanta- 
neous burst is that the bolometric luminosity of the burst 
is now ~ 20% of the total bolometric luminosity of the 
galaxy. 

The mechanical luminosity injected by the SN in the 
"instantaneous" burst (which dominates also in the dou- 
ble burst model) is ~ 10 8 5 L Q over ~ 10 74 yr. This results 
in a total injected energy of ~ 10 erg which is more 
than enough to acco unt fo r the observed superwind. In- 
deed Heckman et al. ( 1990 ) estimate an energy content of 
the winds blown cavity of ~ 1.5 x 10 55 erg (after rescal- 
ing for the different adopted distance of NGC 4945). Both 
models agree with the constraints imposed by dynamical 
measurements that the centr al mas s in stars must be less 
than 6.6xlO 8 M (Koornneef 1993, after rescaling for the 



different assumed distances of the galaxy): the continuous 
SFR would require 5 x 10 9 yr to produce that mass of stars. 

In conclusion, two different star formation histories can 
reproduce the observed starburst properties but only in 
one case does the starburst dominate the bolometric lu- 
minosity of the galaxy. Therefore the available data do not 
allow any constraints on the bolometric luminosity of the 
starburst. 

As shown in the next section, the observed (Lfir/^x) 
ratio of NGC 4945 is equal to that of a "normal" AGN in 
which the £fir is reprocessed UV radiation. If the Lfir in 
NGC 4945 is actually dominated by the starburst, there- 
fore, it is clear that the AGN must be strongly deficient 
in UV relative to X-rays. 

In this starburst-dominant scenario for NGC4945, with 
the black hole mass inferred from the H2O maser mea- 



surements (Greenhill. ct al. 1997), the AGN is emitting at 

< , 



40% of its Eddington Luminosity. 



4-2. NGC 4945 as an AGN dominated object 

By fitting the simultaneous 0.1-200 keV spectrum from 
BeppoSAX, the absorption corrected luminosity in the 
2-10 keV band i s L x (2 - lOkeV) = 3 x lO^ergs" 1 
(Guainazzi et al., |2000| ). If the AGN in NGC 4945 has an 
intrinsic spectral energy distribution similar to a quasar, 
then L x (2-10keV)/L b ol - 0.03 (Elvis et al. |l994j ) there- 
fore (Lboi)AGN ~ 10 44 ergs _1 = 2.6 x 10 10 L m which is the 



10 44 ergs 1 

total far-IR lumin osity of NGC 4945, measured by IRAS 
(Rice et al., |1988|) . Thus, a "normal" AGN in NGC 4945 
could in principle power the total bolometric luminosity. 

For this scenario, we compare NGC 4945 with a nearby 
obscured object, the Circinus galaxy, now considered an 
exam ple of a "sta ndard " Seyfert 2 gala xy (c.f . Oliva et al. 
1994 Oliva et al. |l99S|, Maiolino et al. |l998a|, Matt et al. 
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Fig. 5. Spectral energy distributions of NGC 4945 (upper 
panel) and Circinus (lower panel). "Stars" are the IRAS 
photometric points (except for the points with the highest 
wavelength which are from baloon-borne observations). 
"K", "L", "M" and "N" are the points in the standard 
photometric bands. The hatched areas labeled as "Star- 
burst" and " AGN" represent the continuum levels derived 
from the ionizing photon rates (see text) emitted by star- 
burst and AGN, respectively. The "lOOkeV" points are 
from X-ray observations. The IR spectrum of Circinus 
(solid line in lower panel) is plotted in the upper panel 
(dotted line) after rescaling to match the 100 keV points. 
The solid line in the upper panel is the same spectrum 
after extinction by an extra Ay = 150 mag (see text for 
details). 



1999, Storchi-Bergmann et al. 199S , Curran et al 



In particular, Oliva, Marcon i fc M oorwood ( 1999a ) and, 



1999|) 



previously, Moorwood et al. ( 1996c ) showed that the total 
energy output from the AGN required to explain the ob- 
served emission line spectrum is comparable to the total 
FIR luminosity, concluding that any starburst contribu- 
tion to the bolometric luminosity is small (~10%). The 
choice of the Circinus galaxy is motivated by the simi- 
lar distance (D=4Mpc), FIR and hard X-ray luminosi- 
ties as NGC 4945 (Lfir ~ 1.2 x 10 10 L o ; Siebenmorgen 
et al., |1997| - L x (2 - lOkeV) - 3.4 - 17 x lO^ergs" 1 ; 
Matt et al., |1999[ ). Note that its L x (2 - 10keV)/LFiR ra- 
tio (~ 0.01 — 0.05) is consistent with the average value for 
quasars (Elvis et al. 



1994) 



The overall spectral energy distributions of NGC 4945 
and Circinus are compared in Fig. |^. The "stars" repre- 
sent the IRAS photometric points (except for the points 
with the largest wavelength which are the 150 fim mea- 
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surements by Ghosh et al. |1992| ). In NGC 4945 the points 
labeled with "K", "L" and "N" are the upper limits de- 
rived from our observations, while in Circinus they repre- 
sent emission from the unresolved nuclear source corrected 



for stellar emission (Maiolino et al. 1998a). The points la- 
beled "100 keV" are from Done et al. |1996| (NGC 4945) 
and Matt et al. 1999 (Circinus). The bars between 13.6 
and 54.4 eV are at a level given by vh v ~ Q(H) < hv >, 
where Q(H) is the rate of H-ionizing photons and < hv > 
is the mean photon energy of the ionizing spectrum. For 
NGC 4945, Q(H) is derived from H recombination lines 
and thus represent the energy which is radiated by the 
young starburst; we assumed < hv >= 16 ev. For Circi- 
nus, the point labeled with "Starburst" is similarly derived 
fromBry emission associated with the starburst (Oliva et 
al. |1994[) while that labeled "AGN" is from the estimate 



made by Oliva, Marconi & Moorwood ( 1999a ). 

In the lower panel, we represent the IR spectrum of 
Circinus by connecting the photometric points just de- 
scribed. We plot this same spectrum as a dotted line in 
the upper panel, rescaling to match the 100 keV points. 
NGC 4945 and Circinus have similar X/FIR ratios: 
vL v (100keV)/L FIR ~ 2 x 10~ 3 for Circinus and ~ 3 x 
10~ 3 for NGC 4945. Note that, at each wavelength, both 
Circinus and NGC 4945 were observed with comparable 
resolution. 

If the AGN in NGC 4945 dominates the luminosity and 
its intrinsic spectrum is similar to that of Circinus, then 
the lack of AGN detections in the near-IR and mid-IR re- 
quire larger obscuration. In particular the non-detection 
of a K band point source or of dilution of the CO fea- 
tures can be used to estimate the extinction at 2.2 /xm. In 
Circinus from Maiolino et al. 1998a] (K band) and Matt 



et al. |1999j we can derive vL v {K)/vL v {lWkeV) ~ 0.5. If 
NGC 4945 has a similar near-IR over hard- X-rays ratio 
then, given vL u (K) / vL u (100keV) < 8 x 10~ 4 (K band 
upper limit from CO and X-ray data from Done et al. 



1996 ), the extinction toward the nucleus is AAk > 7mag 
(i.e. AAy > 70mag) larger than in the case of Circi- 
nus. Hot dust in NGC 4945 must be hidden by at least 
Ay > 135 mag, in ag reement with the estimate by Moor- 
wood & Glass ( 1984 ), Ay > 70 and, more recently, with 
an analysis of ISO CVF spectra implying Ay ~ 100 mag 
(Maiolino et al., 2000, in preparation). We note that the 
required extinction is not unexpected and in agreement 
with the X-ray measurements. The measured column den- 
sity in absorption in the X-rays is Nr ~ few x 10 24 cm~ 2 
therefore the expected Ay, assuming a galactic gas-to- 
dust ratio is: 



A 



v 



450 



10 24 cm" 2 



(3) 



The Ay measured from optical/IR data is estimated 
smaller than derived from X-rays (Ay (IR) ~ 0.1 — 



0.5 Ay(X); Granato, Danese & Franceschini 1997 ), there- 
fore the X-ray absorbing column density is in excellent 



agreement with the required extinction. Very high extinc- 
tion, expected in the frame of the unified AGN model are 
observed in many objects as discussed and summarized, 
for instance, in Maiolino et al. 
^999 . 



1998b and in Risaliti et al. 



The higher extinction can also qualitatively explain the 
redder colors of the NGC 4945 FIR spectrum. The solid 
line in the upper panel is the spectrum of Circinus af- 
ter applying foreground extinction by Ay — 150 mag. We 
have applied the extinction law by Draine & Lee ( 1984 ) 
and the energy lost in the mid-IR has been reprocessed as 
40 K dust emission (i.e. black body emission at 40 K cor- 
rected for A -1 ' 75 emissivity). Though a careful treatment 
requires a full radiation transfer calculation, this simple 
plot demonstrates that (i) the redder color of NGC 4945 
with respect to Circinus can be explained with extra ab- 
sorption and (ii) that this is not energetically incompatible 
with the observed FIR luminosity, i.e. the absorbed mid- 
IR emission re-radiated in the FIR does not exceed the 
observed points. 

If the FIR emission is powered by the AGN this is UV 
radiation re-processed by dust. However, if the AGN emits 
~2x 1O 1O L0 in UV photons, high excitation gas emission 
lines should also be observed. The absence of high ioniza- 



tion lines like [O III] A5007, 4959 A (Moo rwood et al. [!996aD 
or [Ne v]A14.3/xm (Genzel et al. 1998| ) and the low exci- 
tation observed in the wind-blown cone strongly argues 
that no ionizing UV photons (i.e. 13.6 < hv < 500 ev) 
escape from the inner region. The low excitation H2/Paa 
map, associated with the peak in H2 emission close to the 
nucleus location, indicates that ALL ultraviolet photons 
must be absorbed within R < l'/5, i.e. R < 30 pc along 
ALL lines of sight. This is in contrast with the standard 
unified model of AGN where ionizing radiation escapes 
along directions close to the torus axis. 

If the AGN is embedded in a thick dusty medium 
then two effects will contribute to its obscuration. First, 
dust will compete with the gas in absorbing UV pho- 
tons which will be directly converted into infrared radi- 



ation (e.g. Netzer & Laor 1993, Oliva, Marconi & Moor- 
wood 1999a). Second, emission lines originating in this 
medium will be suppressed by dust absorption. To es- 
timate the amount of required extinction, note that in 
Circinus [Nev]14.3/rai/[Neii]12.8/xm = 0.4 (extinction 
corrected) and in NGC 4945 [N e v] /[ Ne ii] < 0.008 (both 



ratios are from Genzel et al. [L998D . If NGC4945 has 
the same intrinsic ratio as Circinus, then the observed 
[Nev]/[Neii] ratio requires A(14.3/im) > 4.2mag corre- 
sponding to Ay > llOmag and in agreement with the 
above estimates. 

We conclude that the AGN can power the FIR emis- 
sion if it is properly obscured. Inferring the black hole mass 
from the H2O maser observations (1.4 x 10 6 M Q , Green- 
hill, et al. 1997 ), we find in this scenario that the AGN is 
emitting at ~ 50% of its Eddington Luminosity. 
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4-3. On the existence of completely hidden Active 
Galactic Nuclei 

As discussed above, if an AGN powers the FIR emission of 
NGC 4945, it must be hidden up to mid-IR wavelengths 
and does not fit in the standard unified model. The pos- 
sible existence of such a class of Active Nuclei, detectable 
only at > 10 keV, would have important consequences on 
the interpretation of IR luminous objects whose power 
source is still debated. 

Genzel et al. fll998| ) and Lutz et al. ( |1998| ) com- 
pared mid-IR spectra of Ultra Luminous IRAS galaxies 
(ULIRGs, see Sanders & Mirabel 1996| for a review) with 
those of AGN and starburst templates. They concluded 
that the absence of high excitation lines (e.g. [Ne v]) and 
the presence of PAH features undiluted by strong thermal 
continuum in ULIRGs spectra strongly suggest that the 
starburst component is dominant. They also show that, 
after a proper extinction correction, the observed star for- 
mation activity can power FIR emission. In their papers, 
NGC 4945 is classified as a starburst because of its mid- 
IR properties but, as shown in the previous section, NGC 
4945 could also be powered by a highly obscured AGN and 
the same scenario could in principle apply to all ULIRGs. 
Their bolometric emission can be powered by an active 
nucleus completely obscured even at mid-IR wavelengths. 

The same argument could be used for the sources de- 
tected at submillimeter wavelengths by SCUBA which can 
be considered as the high redshift counterpart of local 
ULIRGs. If they are powered by hidden active nuclei then 
their enormous FIR emission would not require star for- 
mation rates in excess of > 100M Q yr _1 (e.g. Hughes et 
al. 1998), and this would have important consequences for 



black holes and the high redshift AGN component (see 



understanding the history of star formation in high red- 
shift galaxies. 

In addition, it is well known that in order to explain 
the X-ray background a large fraction of obscured AGN 
is required. However Gilli et al. ( |1999 ) have shown that, 
in order to reconcile the observed X-ray background with 
hard X-ray counts, a rapidly evolving population of hard 
X-ray sources is required up to redshift ~ 1.5. No such 
population is known at the moment and the only class 
of objects which are known to undergo such a ra pid den- 



sity evolution are local ULIRGs (Kim et al. 1998 ) and, a t 



higher redshift, the SCUBA sources (Smail et al. 1997 ) 
SCUBA sources are therefore candidates to host a popu- 
lation of highly obscured AGNs. 



Almaini et al. (1999) suggest that, if the SED of high 
redshift AGN is similar to those observed locally, one 
can explain 10-20% of the 850 /im SCUBA sources at 1 
mJy. This fraction could be significantly higher if a large 
population of AGN are Compton thick at X-ray wave- 
lengths. Trentham, Blain & Goldader (1999) show that 



if the SCUBA sources are completely powered by a dust 
enshrouded AGN then they may help in explaining the 
discrepancy between the local density in super massive 



also Fabian & Iwasawa 1999). 

Establishing the nature of SCUBA sources could be ex- 
tremely difficult if the embedded AGNs are like NGC4945, 
i.e. completely obscured in all directions, because they 
would then not be identifiable with the standard opti- 
cal/IR diagnostics. Incidentally, this fact could possibly 
account for the sparse de tectio ns of type 2 AGNs at high 
redshifts (Akiyama et al. 1999 ). 

The best possibility for the detection of NGC4945- 
like AGNs is via their hard X-ray emission but, unfortu- 
nately, the sensitivity of existing X-ray surveys is still not 
high enough to detect high z AGN and the low spatial 
resolution makes identifications uncertain in the case of 
faint optical/near-IR counterparts. Moreover, hard X-rays 
alone are not enough to establish if the AGN dominates 
the bolometric emission. 



5. Conclusions 

Our new HST NICMOS observations of NGC 4945, com- 
plemented by new ground based near and mid-IR observa- 
tions, have provided detailed morphology of the nuclear re- 
gion. In Paa, we detect a lOOpc-scale starburst ring while 
in H2 we trace the walls of a conical cavity blown by su- 
pernova driven winds. The continuum images are strongly 
affected by dust extinction but show that even at HST res- 
olution and sensitivity, the nucleus is completely obscured 
by a dust lane with an elongated, disk-like morphology. 
We detect neither a strong point source nor dilution in 
CO stellar features, expected signs of AGN activity. 

Whereas all the infrared properties of NGC 4945 are 
consistent with starburst activity, its strong and variable 
hard X-ray emission cannot be plausibly accounted for 
without the presence also of an AGN. Although the star- 
burst must contribute to the total bolometric luminosity 
we have shown, using starburst models, that the actual 
amount is dependent on the star formation history. A 
major contribution from the AGN is thus not excluded. 
Irrespective of the assumption made, however, our most 
important conclusion is that the observed variable hard 
X-ray emission combined with the lack of evidence for 
reprocessed UV radiation in the infrared is incompatible 
with the "standard" AGN model. If the AGN dominates 
the bolometric luminosity, then its UV radiation must be 
totally obscured along all lines of sight. If the starburst 
dominates then the AGN must be highly deficient in its 
UV relative to X-ray emission. The former case clearly 
raises the possibility that a larger fraction of ULIRGs than 
currently thought could actually be AGN rather than star- 
burst powered. 

Acknowledgements. A.M. and A.R. acknowledge the partial 
support of the Italian Space Agency (ASI) through grants 
ARS-98-116/22 and ARS-99-44 and of the Italian Ministry for 
University and Research (MURST) under grant Cofin98-02-32. 



14 



Marconi et al.: The elusive AGN of NGC 4945 



E.J.S. acknowledge support from STScI GO grant O0113.GO- 
7865. We thank Roeland P. van der Marel for the use of the 
pedestal estimation and quadrant equalization software. 

References 

Akiyama M., et al., 1999, Proceeding of the first XMM wo rk- 
shop "Science with XMM", in press (|astro-ph/9811012|) 



1999, MNRAS, 305, L59 



Almaini O., Lawrence A., Boyle B.J. 
Baan W.A., 1985, Nature, 315, 26 

Biretta J. A., et al., 1996, WFPC2 Instrument Handbook, Ver- 
sion 4.0 (Baltimore:STScI) 

Blain A.W., Kneib J.P., Ivison R.J., Smail I., 1999, ApJ, 512, 
L87 

Brock D., 1988, ApJ, 329, 208 

Bushouse H., Skinner C.J., MacKenty J.W., 1997, NICMOS 

Instrument Set. Kept. 97-28 (BaltimoreiSTScI) 
Cappi M., et al., 1999, A&A, 350, 777 

Comastri A., Setti G., Zamorani G., Hasinger G., 1995, A&A, 
296, 1 

Curran S.J., et al., 1999, A&A, 344, 767 
Done C, Madejski G.M., Smith D.A., 1996, ApJ, 463, L63 
Dos Santos P.M., Lepine J.R.D., 1979, Nature, 278, 34 
Draine B.T., Lee H.M., 1984 ApJ, 285, 89 
Elmouttie M., et al., 1997, MNRAS, 284, 830 
Elvis M., et al., 1994, ApJS, 95, 1 
Fabian et al., 1998, MNRAS, 297, Lll 
Fabian A.C., Iwasawa K., 1999, MNRAS, 303, L34 
Forbes D.A, Norris R.P, 1998, MNRAS, 300, 757 
Genzel R., Lutz D., Sturm E., Egami E., Kunze D., Moor- 
wood A.F.M., Rigopoulou D., Spoon H.W.W., Sternberg 
A., Tacconi-Garman L.E., Tacconi L., Thatte N., 1998, 
ApJ, 498, 579 
Ghosh S.K., et al., 1992, ApJ, 391, 111 
Gilli R., Risaliti R., Salvati M., 1999, A&A, 347, 424 
Granato G.L., Danese L., Franceschini A., 1997, ApJ, 486, 14 
Greenhill L.J., Moran J.M., Herrnstein J.R., 1997, ApJ, 481, 
L23 

Guainazzi M., et al, 2002000, A&A 



ph/3001528) 



in press ( astro 



1997 ApJS, 108, 



P., Mestreau P., 
35, 203 



Heckman T.M., Armus L., Miley G.K., 1990, ApJS, 74, 833 
Hughes D.H., et al., 1998, Nature, 394, 241 
Hunt L.K., Malkan M.A., Salvati M., et al.. 
229 

Iwasawa K., et al., 1993, ApJ, 409, 155 
Kaufl H.U., Jouan R., Lagage P.O., Masse 
Tarrius A., 1994, Infrared Phys. Technol. 
Kim D.-C, Sanders D.B., 1998, ApJS, 119, 41 
Koornneef J., 1993, ApJ, 403, 581 

Koornneef J., Israel F.P., 1996, New Astronomy, 1, 271 
Leitherer C, et al., 1999, ApJS, 123, 3 

Lutz D., Spoon H.W.W., Rigopoulou D., Moorwood A.F.M., 

Genzel R., 1998, ApJ, 505, L103 
MacKenty J.W., et al., 1997, NICMOS Instrument Handbook, 

Version 2.0 (Baltimore:STScI) 
Maiolino R., Ruiz M., Rieke G.H., Keller L.D., 1995, ApJ, 446, 

561 

Krabbe A., Thatte N., Genzel R., 1998, ApJ, 493, 



Maiolino R. 

650 
Maiolino R 

Matt G 



Maiolino R., et al., 1999, ApJ, in press ( |astro-ph/9910160| ) 
Marconi A., Schreier E.J., Koekemoer A., Capetti A., Axon 

D.J., Macchetto F.D., Caon N., 2000, ApJ, 528, 276 
Matt G., Guainazzi M., Maiolino R., Molendi S., Perola G.C., 
Antonelli L.A., Bassani L., Brandt W.N., Fabian A.C., 
Fiore F., Iwasawa K., Malaguti G., Marconi A., Poutanen 
J., 1999, A&A, 341, L39 
Mauersberger R., et al., 1996, A&A, 309, 705 
Moorwood A.F.M. and Glass I.S., 1984, A&A, 135, 281 
Moorwood A.F.M. and Oliva E., 1988, A&A, 203, 278 
Moorwood A.F.M., Finger G., Biereichel P., Delabre B., Van 
Dijsseldonk A., Huster G., Lizon J. L., Meyer M., Gemper- 
lein H., Moneti A., 1992, The Messenger, 69, 61 
Moorwood A.F.M. and Oliva E., 1994a, ApJ, 429, 602 
Moorwood A.F.M., Finger G., Gemperlein H., 1994b, The Mes- 
senger, 77, 8 

Moorwood A.F.M., van der Werf P.P., Kotilainen J.K., Mar- 
coni A., Oliva E., 1996a, A&A, 308, LI 
Moorwood A.F.M., Space Sci. Rev., 1996b, 77, 303 
Moorwood A.F.M., Lutz D., Oliva E., Marconi A., Netzer H., 
Genzel R., Sturm E., de Graauw Th., 1996c, A&A, 315, 
L109 

Netzer H., Laor A., 1993, ApJ, 404, L51 

Oliva E., Salvati M., Moorwood A.F.M., Marconi A., 1994, 
A&A, 288, 457 

Oliva E., Origlia L., Kotilainen J.K., Moorwood A.F.M., 1995, 
A&A, 301, 55 

Oliva E., Marconi A., Cimatti A., di Serego Alighieri S., 1998, 
329, L21 

Oliva E., Marconi A., Moorwood A.F.M., 1999a, A&A, 342, 87 
Oliva E., Origlia L., Maiolino R., Moorwood A.F.M., 1999b, 

A&A, 350, 9 
Persic M., et al., 1998, A&A, 339, L33 
Rice G.H., et al., 1988, ApJS, 68, 91 
Richstone D., et al., 1998, Nature, 395, 14 
Rieke G.H., et al., 1980, ApJ, 238, 24 
Risaliti G., et al., 1999, ApJ 522, 157 
Rowan-Robinson M., et al., 1997, MNRAS, 298, 490 
Sanders D.B., Mirabel I.F., 1996, ARA&A, 34, 749 
Schreier E.J., Marconi A., Axon D.J., Caon N., Macchetto D., 

Capetti A., Hough J.H., Young S., Packham C, 1998, ApJ, 

499, L143 (Paper II) 
Siebenmorgen R., Moorwood A., Freudling W., Kaeufl H.U., 

1997, A&A, 325, 450 
Skinner C.J., Bergeron L.E., Daou D., 1998, HST Calibration 

Workshop, ed. S. Casertano et al. (Baltimore:STScI), in 

press 

Smail I., Ivison R.J., Blain A.W., 1997, ApJ, 490, L5 
Storchi-Bergmann T., et al., 1999, MNRAS, 304, 35 
Trentham N., Blain A.W., Goldader J., 1999, MNRAS, 305, 61 
Voit, M., et al., 1997, HST Data Handbook - Vol. I, Version 3.0 
(BaltimoreiSTScI) 



Salvati M., Bassani L., Dadina M., della Ceca R.. 
Risaliti G., Zamorani G., 1998, A&A 338, 781 



